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ABSTRACT 

Aims. This is the third paper of a series devoted to study the properties of bars from long slit spectroscopy to understand their 
formation, evolution and their influence on the evolution of disk galaxies. In this work we aim to determine the gas metallicity 
distribution of a sample of 20 barred early-type galaxies. We compare the nebular and stellar metallicity distributions to conclude 
about the origin of the warm gas. 

Methods. Long-slit spectroscopy along the bar was obtained and metallicities derived using different calibrations. We compare the 
results of nebular emission metallicities using different semi-empirical methods. We carry out AGN diagnostic diagrams along the 
radius to determine the radius of influence of the AGN and the nuclei nature of the studied galaxies. We then derive the gas metallicities 
along the bars and compare the results to the distribution of stellar metallicities in the same regions. 

Results. Most of the gas emission is centrally concentrated, although 15 galaxies also show emission along the bar. In the central 
regions, gas oxygen abundances are in the range 12+log(0/H)= 8.4-9.1. The nebular metallicity gradients are very shallow in the 
bulge and bar regions. For three galaxies (one of them a LINER), the gas metallicities lie well below the stellar ones in the bulge 
region. These results do not depend on the choice of the semi-empirical calibration used to calculate the abundances. We see that the 
galaxies with the lowest abundances are those with the largest rotational velocities. Unlike with the stellar abundances, we find no 
correlation between the nebular abundances and the central velocity dispersion. In most galaxies the slope for the nebular abundance 
distribution in the bulge region is shallower than that for the stellar metallicity. 

Conclusions. The presence of gas of significantly lower metallicity than the stellar abundances in three of our galaxies, points to an 
external origin as the source of the gas that fuels the present star formation in the centre of some early-type barred galaxies. The fact 
that the bar/disk nebular metallicities are higher than the central ones might be indicating that the gas could be accreted via cooling 
flows instead of radial accretion from gas sitting in the outer parts of the disk. 
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^SJ ' 1. Introduction tories compared to the bulges of their unbarred counterparts (see 

\Q , also Ellison et al. 2001 for a similar conclusion for gas-phase 

■ With the exception of few primordial light elements, stellar nu- metallicities). 
^ ; cleosynthesis is responsible of the secular metal enrichment Gas and stars respond to the gravitational potential due to 

^ , in galaxies. The final metallicity distribution within a galaxy the i r OW n nature, being viscosity and magnetic fields only im- 

t-H ' is reasonably well reproduced by chemical evolution models por tant in the former component. The highly asymmetric bar 

; considering the appropiate star formation history together with potential induces differential radial motions. While the gas com- 

. ^ , episodes of gaseous inflow and outflows (e.g. Portinari & Chiosi pone nt, being highly dissipative, suffers from the gravitational 

K> ; 1999; Chiappini et al. 2000; Prantzos 2008, Colavitti et al. 2009). tor q Ue of the non-axisymmetric mass component, the stars are 

^ i However, the metallicity gradients can be futher modified by dy- ma inly affected by orbital mixing. Furthermore, studies of the 

_ C3_. namical processes (e.g., Roskar et al. 2008, Schonric & Binney gas -phase abundances provide with present-day snap-shots of 

2009; Sanchez-Blazquez et al. 2009). It is, therefore, evident that me interstellar medium abundance. On the other hand, the study 

observing how metals are distributed in a galaxy should highly G f ^Uar abundances provide archeological clues as to the for- 

constrain its evolution. However, the detailed dynamical pro- mat i on an d evolution of the bar. Therefore, the study of stellar 

cesses responsible for the modification of the metallicity distri- an d gas metallicities are of crucial importance to interpret the 

bution and their importance are not yet well understood. processes dominating galaxy evolution. It is expected that, if gas 

In particular, it is not known, from the observational point is injected into the interstellar medium (ISM) from stellar out- 

of view, the importance that bars might have in producing the puts, it would generally be more metal-rich than if it has an ex- 

metallicity radial mixing. Bars are believed to affect the overall ternal origin. 

dynamics of the galaxy and are a well known mechanism to in- Deriving the disk radial distribution of stellar abundances 

duce secular evolution (Athanassoula 2003; Pfenniger & Friedli from spectroscopic measurements is difficult due to the low sur- 

1991). For instance, in a previous work (Perez & Sanchez- face brightness of the disk and the contamination of the emis- 

Blazquez 201 1, hereinafter, Paper II) we hinted that the bulges of sion lines from ionised gas. Few works have tried to overcome 

early-type barred galaxies show different stellar enrichment his- these difficulties and have derived the radial distribution for a 
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small number of galaxies (e.g. Yoachim & Dalcanton 2008; 
MacArthur, Gonzalez & Courteau 2009; Sanchez-Blazquez et 
al. 2011). The stellar metallicity gradients along the bars of the 
20 early-type galaxies presented here were carried out in previ- 
ous works (Perez et al. 2009, hereinafter, Paper I) and provided 
interesting results regarding the radial distribution of the stellar 
parameters in the bar region. 

As for the nebular gas abundance distribution, for most spi- 
ral galaxies, both barred and unbarred, negative radial metal- 
licity gradients along the disk have been obtained for the gas 
component (e.g. Bresolin et al. 2009). These gradients seem to 
be shallower in the case of barred galaxies than in unbarred 
galaxies, and larger for late-types than for early ones (Pag el & 
Edmunds 1981; Alloin et al. 1981; Vila-Costas and Edmunds 
1992; Zaritsky et al. 1994). The issue is not closed yet, so work is 
being done to investigate the variation of the gradient with radius 
(Bresolin et al. 2009 obtain a higher gradient in the central part 
of M83, and flatter in the outer disk, beyond 1.2 R25) and also 
in azimuth (Balser et al. 201 1, for the Milky Way). Simulations 
explain the shallower gradient value in the outer part than in the 
inner side of a galaxy based on radial mixing processes associ- 
ated to a strong bar (e.g. Considere et al. 2000; Zahid & Bresolin 
201 1; Friedli et al. 1994; Friedli 1999). In this work, we will fo- 
cus on the analysis of the gas metallicity, in the bar region of the 
galaxies presented in Papers I and II. 

Only a few works have considered both, nebular and stellar 
metallicities. A pioneering paper in this study is that of Storchi- 
Bergmann et al. (1994). They found a good correlation between 
stellar and gas metallicity investigating through the correlation 
between the oxygen abundance in the nebular component, and 
the absorption-line equivalent width W(CIV/11550) (as a tracer 
of stellar metallicity) in a sample of 44 star-forming galaxies. 
They concluded that galaxies with lower values of metallicity 
appear to experience instantaneous outbreaks of star formation, 
but those with higher metallicity are probably experimenting an 
ongoing stellar formation. Annibali et al. (2010) estimated gas 
abundance gradients in early-type galaxies, and they compared 
these results with their stellar metallicities (Annibali et al. 2007), 
concluding that the gas metallicity tends to be lower than the 
stellar one, beeing this effect higher for higher metallicities. 

In this paper, we present an analysis of the radial oxygen 
abundance distribution in the bulge and bar region of a sample of 
20 early type barred galaxies, with and without nuclear activity. 
The stellar component of this sample has been studied in Papers 
I and II. In Sec. 2 observations and data reduction are described. 
The method to analyse the spectral lines to estimate their flux is 
presented in Sec. 3. In Sec. 4 we show the diagnostic diagrams 
and derive the abundance distributions with different methods. 
A comparison of the nebular and stellar metallicity gradients is 
also presented in this section. The results are discussed in Sec. 
5. 

2. Observations and data reduction 

We obtained long-slit spectra along the bar major axisS for 
our sample of 20 barred galaxies. The observations were per- 
formed in two different runs, with the double beam spectro- 
graph at Siding Spring Observatory (Australia, hereinafter, runl) 
and with the IDS spectrograph at the Isaac Newton Telescope 
(La Palma, Spain, hereinafter, run2). The observations were de- 
scribed in detail in Paper I. To summarise, spectra for seven 

1 The bar position angles were derived using the Digital Sky Survey 
(DSS) images. 



galaxies were obtained in the first run, covering a wavelength 
range from 3892-5815A and 5390-7314A with a spectral reso- 
lution of FHWM~2.2A. In spite of this overlap between both 
ranges, relative calibration is not possible because the trans- 
mission of the dichroic gets as low as 50%, decreasing the 
S/N relation, and we have no bright lines in this region. In the 
second run, spectra for 13 galaxies were obtained, covering a 
wavelength range of 3020-6665A with a spectral resolution of 
~3A (FWHM), but the unvignetted range is 3600-6300A. 

The seeing size during these observations range from 0.6" to 
1.2" at the INT and from 1" to 1.5" at SSO. The characteristics 
of the sample are described in Paper I and II. For convenience, 
Table[T]shows the sample main properties. 

Atmospheric dispersion can have important effects in spec- 
trophotometry when the slit position angle differs from the paral- 
lactic angle during the observations. This is not so critical in our 
observations as we observed extended objects (diffuse emission 
or emission knots with sizes typically a few arcseconds larger 
than our slit width). 

In any case, to be in the safe side, we kept the target airmass 
below 1.6 during the observations for ensuring an atmospheric 
differential refraction between [OIIJ3727A and [OIII]5007A 
smaller than 1.5" (Filippenko 1982). Differential atmospheric re- 
fraction between any other pair or lines involved in the empirical 
metallicity calibrations used in the paper are always smaller than 
the differential atmospheric refraction between [Oil] 3727 A and 
[OIII]5007A. 

For a few galaxies in which the difference between the bar 
position angle and the parallactic angle was greater than « 40 
deg, we tried to keep the airmass below 1.3 (which implies 
a differential atmospheric refraction between [Oil] 3727 A and 
[OIII]5007A of « 1"). 

The reduction of the two runs was carried out with the pack- 
age REDUCEME (Cardiel 1999) and IRAF tasks. Standard data 
reduction procedures (flat-fielding, cosmic ray removal, wave- 
length calibration and sky subtraction) were performed. Error 
images were created at the beginning of the reduction and were 
processed in parallel with the science images. For details about 
the reduction steps, see Paper I. 

3. Analysis 

We have observed 20 early type barred galaxies. Nebular emis- 
sion has been detected in 19 of them. In NGC 2950, we detect 
no emission at all, and in NGC 4245 only at a galactocentric 
radius of 6 arcsec. In NGC 1832 and NGC 1530 we resolve 
knots of star formation (SF hereafter) at radius R< 0.25 R25 (R25 
being the length of the projected major axis of a galaxy at the 
isophotal level 25 mag/arcsec 2 in the B-band). For the other 16 
galaxies, we observe non-resolved emission extending at least 
along the bulge radius, i.e. R < 5" - 15" depending on the 
galaxy (Table I). For seven galaxies (NGC 1433, NGC 1530, 
NGC 1832, NGC 2665, NGC 2935, NGC 3081 and NGC 4314) 
we have also found knots of gas emission at larger radii, corre- 
sponding to one or several crowded H 11 regions. 

To compare gas and stellar abundances we extract the spectra 
along the radius using the same binnings than in Paper I. These 
binnings were chosen to ensure errors lower than 15% in most 
of the Lick/IDS indices. We also add apertures at larger radius to 
include those knots of star formation found in 7 of our galaxies 
(see above). 

The gas emission line fluxes are measured from the spec- 
tra after subtracting stellar templates. These stellar templates are 
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Table 1. General properties of the galaxy sample. 
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(1) Galaxy identification name. (2) Galaxy morphological classification. (3) Bar class derived from: " the A* - band light distribution, Block et al. 
(2001); * the K - band light distribution, Block et al. (2004); the K - band light distribution, Buta et al. (2006); d the H - band light distribution, 
Laurikainen et al. (2004). (4) Nuclear type, from this paper when possible. Between "()" is the previous classification carried out by Veron-Cetty & 
Veron (2006). T is transition object. "-" means that no classification as active galaxy has been found in the literature. (5) Inner morphology, from 
" Buta (1986) b Jungwiert et al. (1997) c Erwin (2004) d Wozniak et al. (1995). (6), (7) & (8) Apparent total B-magnitude, maximum rotational 
velocity corrected for inclination and galaxy inclination, from the Hyperleda galaxy catalog. (9) Bar semi-major axis, form Paper II. (10) Bulge 
size, from Paper II. (11) runl are galaxies observed at the SSO and run2 those observed at the INT. 



obtained using GANDALF (Sarzi et al. 2005) and the set of stellar 
population models by Vazdekis et al. (2010) using the MILES li- 
brary (Sanchez-Blazquez et al. 2006). After the stellar template 
subtraction, the spectra contains only the emission from the gas, 
with no contribution from the underlying stellar population. As 
representative examples, we show in Figs. Q] [2] and [3]the central 
spectrum for three galaxies, two of run2 (one of them with weak 
emission and the other with stronger emission lines) and one of 
runl. They are shown before and after the template subtraction. 

The emission line fluxes are measured using the SPL0T task 
in IRAF, by integrating the line intensity over a locally fitted 
continuum level. After marking two continuum points at both 
sides of the line, the flux is obtained taking into account the area, 
without a previous fitting to a Gaussian or to any other assumed 
profile. Therefore asymmetries due to dynamical processes are 
indirectly considered. The corresponding uncertainties are esti- 
mated adding in quadrature statistical errors measured by SPL0T, 
flatfielding errors and the error of the flux calibration. Typical er- 
rors are of order ~5% (Table [2] shows the errors for the central 
apertures). 

The observed line ratios relative to H/3 are afterwards cor- 
rected for interstellar reddening, by determining the reddening 
coefficient c(H/3) through the expression: 

_H> = g(j> 1 QCW)[f W -WP)] , J ) 

1{HP) F(Hp) 

where f(A) is the Seaton (1979) reddening curve parametrized 
by Howarth (1984), and 1(A) and F(A) are the intrinsic and ob- 
served emission line flux at wavelength A, respectively. We em- 
ploy the Balmer lines detected within our spectral ranges. As ex- 
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Fig.l. Central spectrum for NGC 1169 (from run2). The top 
panel shows the observed spectrum (yellow line) and the fitted 
template (black line). Bottom panel shows the residuals from the 
fit. Oxygen and H/3 lines are clearly visible. 



plained in Sect. 2, both arms of the runl observations have been 
treated separately for avoiding the introduction of an additional 
continuum scaling error, and therefore the Ha/H/3 is not avail- 
able. The extinction constant has then been obtained from the 
average of the extinction constant obtained from the fitting of all 
available Balmer line ratios relative to H/3, and the corresponding 
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Fig. 2. Central spectrum for NGC 2273 (run2). Panels are the 
same as in Fig. [T] In this spectrum it can be also identified the 
[Nem]^3868 line. 
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Fig. 3. Blue and red central spectra for NGC 3081 (runl). In the 
top pannels observed spectra (yellow lines) and templates (black 
lines) are shown. Bottom panels show the residuals of the fit. 



intrinsic or theoretical values. We consider the intrinsic flux ratio 
from Hummer & Storey (1987), assuming T e = 10 4 K and N e = 
100 crrT 3 . For those cases in which there is a high dispersion 
between the values for c(HyS) obtained from different lines ra- 
tios, we use the extinction constant obtained from HyfH/3 alone. 
We assign an internal extinction of zero to objects for which the 
observed Balmer decrement is lower than the theoretical values. 
The derived reddening is negligible for most of the spectra. In 
the case of NGC 4245, we have obtained only H/3 and no other 
Balmer line, but we have not considered this galaxy in the dis- 
cussion. For NGC 1 169, we observe only H/3 except at 5 galac- 
tocentric radii. For 4 galaxies the observed Balmer decrement is 
lower than the theoretical one, so we consider as zero the internal 
extinction. In the fifth case, the obtained internal extinction is so 
low that differences between line fluxes considering or not this 
internal extinction are lower than the obtained errors. Therefore, 
we have assumed a negligible extinction in this galaxy. 

Emission line strengths for the central aperture are shown 
in Table [2] For NGC 4314 we show the emission line strengths 



for the innermost aperture (0.8 arcsec from the galaxy center) 
with detected emission. This central aperture corresponds to 0.75 
arcsec for runl galaxies and 0.4 arcsec forrun2 galaxies. Values 
of [On], [Nem] and [Om] are relative to H/3=100, and those of 
[Nil] and [S n] to Ha= 100. 

For those galaxies presenting knots of gas emission, we also 
extract spectra with apertures including all detected emission. In 
these cases, we correct from extinction using an iterative method 
to determine, simultaneously, the underlying stellar absorption 
and the interstellar reddening: we derive the extinction coeffi- 
cient c(H/3) for an assumed equivalent width of the stellar ab- 
sorption EW fl f,j (the same for all observed Balmer lines). The 
assumed EW,,;,, is varied until convergence between the c(H/3) 
obtained from all the Balmer line ratios in use. The c(H/3) values 
are in the range 0-2, and EW„£, was found to be in the range 
0.7-1.8. 



4. Results 

4.1. Star formation versus nuclear activity 

Once the line fluxes are measured and extinction corrected, 
we can represent diagnostic diagrams commonly used to dis- 
tinguish between star formation (SF) regions and regions con- 
taining or affected by an active nucleus (AGN). However, clas- 
sical diagnostic diagrams (log([Oiii]/H/?) versus log([Nii]/Ha') 
and log([Om]/Hj0) versus log([S n]/Ha)) (Baldwin et al. 1981, 
Veilleux & Osterbrock 1987) can only be applied to runl, due 
to the limited spectral coverage of run2. For four galaxies, one 
from runl and three from run2 we detect the [Ne m]/13868 line. 
This line has been proposed as an empirical indicator of metal- 
licity (Nagao et al. 2006) and also as a diagnostic to distinguish 
between starburst galaxies and AGNs (see discussion in Perez- 
Montero et al. 2007). 

Figures [4] [5] and |6]show the diagnostic diagrams. Figures 
and [5] show these diagnostic diagrams at different galactocen- 
tric distances. In Fig. [4] we have plotted the theoretical curves 
of Kewley et al. (2001), which sets the maximun expected 
[O m]/H/3 for a given [N n]/Ha for a star forming region, those 
from Kauffmann et al. (2003) which sets the empirical division 
between star-forming galaxies and transition objects, and that of 
Kewley et al. (2006) separating Seyfert and LINERs galaxies. In 
Fig. [5] we plot the relation between the line ratios log([S ii]/Hq-) 
and log([Ora]/H/3) for the galaxies of our sample at different 
radii. The theoretical models separating the values correspond- 
ing to ionization by star formation and AGN are also indicated. 
Rola et al. (1997) proposed an alternative diagnostic diagram 
based on the [Nem]/13869 line. In Fig.|6]this diagnostic diagram 
for galaxies in which we have detected this line is plotted, re- 
lating log^NemVHyS) with log([Om]/H/3). The solid line indi- 
cates the limit between star formation regions and AGNs, and the 
dashed one between LINERs ([Om]5007/Hj8< 0.5) and Seyferts 
2 ([Oin]5007/HyS> 0.5) as explained in Rola et al. (1997). 

Figures|4]and [5]also show that at large radii there is a larger 
probability that the gas is ionized by star forming regions, while 
central emission is characterized by AGN-like line-ratios. 

Studying the diagnostic diagrams for individual galaxies, 
we confirm the classification carried out by Veron-Cetty & 
Veron (2006) with the exception of NGC 5101 and NGC 1433. 
According to the emission line in the central region, we reclas- 
sify NGC 5101 as a LINER even if it was not classified in 
such a way in Veron-Cetty & Veron (2006), but in agreement 
with Moiseev (2001). The runl galaxies that are not classified 
as Seyfert or LINER, NGC 2665 and NGC 2935, are transition 
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Table 2. Emission Line Strengths for the central aperture. Emission line strengths for the central aperture of the galaxies, except for 
NGC 43 14, which is at 0.8 arcsec from the galaxy center. The line fluxes of [O n], [Ne m] and [O in] lines are relative to Hf3 of 100 
and those of [N n] and [S n] lines are relative to Ha of 100. 
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Fig. 4. Diagnostic diagram for the emission-line regions of runl 
galaxies at different galactocentric distances (R) scaled with R25. 
The solid line (Kewley et al. 2001) sets the star formation up- 
per limit, the dashed line (Kauffmann et al. 2003) the empirical 
division between star-forming galaxies (SF) and transition ob- 
jects (T), and the horizontal line (Kewley et al. 2006) separates 
Seyferts and LINERs. 
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Fig. 5. Diagnostic diagram for the emission-line regions of runl 
galaxies for different galactocentric distances (R) scaled with 
R25. The solid line establishes the separation between star forma- 
tion regions and AGNs, below the dashed one there are LINERs 
and above it Seyferts (Kewley et al. 2006). 



objects. In Table I we show the nuclear type obtained in this pa- 
per and in previous classification carried out by Veron-Cetty & 
Veron (2006). 

Since we have radial information for the diagnostic dia- 
grams, we have considered the question of the AGN influence 
region by analyzing the radial behavior of the diagnostic ratios, 
to find the radius at which star formation starts being the main 
photoionisation source. Fig. [7] shows the radial distribution of 
log([Om]/HjS) (or log([Nein]/Hy6)). In the same figure we have 
considered the measured values of [N ii]/Hq'. With these values 



and taking into account the theoretical curve of Kewley et al. 
(2001) (or Rola et al. 1997 for log([Ne iii]/H/?)) we obtain the 
limiting values to consider star formation or nuclear activity. 
They are indicated with a dashed line. At a given galactocentric 
radius, if the observational points are above it, a star forming re- 
gion is not the main photoionization component at that radius. 
There are differences among galaxies. The AGN or LINER is 
the dominant ionization mechanism up to different radii for the 
different galaxies: For example, for NGC 2273 it is just AGN 
ionized at the lowest radius point. For NGC 2217 the influence 
region of AGN extends to ~ 0. 1 kpc. 
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Fig. 7. Radial distribution of log([Ora]/Hj8) (for active galaxies of runl) and log([Ne m]/H/J) (forrun2 galaxies with [Nem]/13869 
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Fig. 6. Diagnostic diagram for those galaxies in which we 
have detected the [Nem]/13869 line. The solid line indicates 
the limit between star formation regions and AGNs, and the 
dashed one between LINERs ([O iii]5OO7/Hj0< 0.5) and Seyferts 
2 ([Om]5007/H/J> 0.5) (Rola et al. 1997). 

4.2. Oxygen abundance and gradients 

The measurement of nebular gas abundances presents sev- 
eral difficulties. The direct method to calculate oxygen abun- 



dances implies the use of ratios between temperature sensi- 
tive lines, such as [O m]H4959,5007M4363 (Osterbrock 1989). 
Unfortunately, when the metallicity increases, the electronic 
temperature decreases and the auroral lines (e.g. [O in]/l4363) 
become too weak to be measured in most extragalactic objects 
and, therefore, we have to find empirical, semiempirical and the- 
oretical calibrations to calculate abundances. Some relations be- 
tween the gas metallicity and line flux ratios of strong emission 
lines have been found. A drawback of these indicators is that 
they have been basically developed for low-metallicity H n re- 
gions, where the [O m]/l4363 line is easily observed. However, a 
considerable effort to study high metallicity H 11 regions has also 
been made (e.g. Castellanos, Diaz & Terlevich 2002, Bresolin 
2007), since these are the most abundant in early spiral galaxies 
and in the central regions of late type spiral galaxies (Vila-Costas 
& Edmunds 1992). Although semiempirical methods provide 
better results in statistical studies than individual values (Perez - 
Montero et al. 2007), we use these methods to estimate the oxy- 
gen abundance variations with galactocentric radius. Systematic 
discrepancies among different methods have been found (see, 
for instance, Kewley & Ellison 2008; Lopez-Sanchez & Esteban 
2010, and references therein), so we use several methods to study 
the dependence of the results on the method chosen. 

An additional problem is that 1 1 out of our 20 galaxies have 
been classified as LINERs or Seyferts. The mechanism of ion- 
ization in the case of LINERs is still not well understood, so 
there are no methods to estimate abundances in the central parts 
of these galaxies. Some methods exist in the case of Seyferts, but 
due to spectral coverage we can only apply one of them, and in 
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some particular cases. We also apply this method in LINERs, 
which may be intermediate between AGN excitation or pho- 
toionization due to star formation. 

4.2.1 . Nebular abundances: Errors in their determination and 
comparison between different methods 

We use different methods to estimate the nebular abundances 
for the two different runs, as the spectral range covered in 
both data sets is different. The brightest lines measured, apart 
from the Balmer lines (Hy, H/3, Ha) are [Om]/U4959, 5007, 
[Nn]^6548, 6583, [Sn]^6717, 6731 for runl galaxies, and 
[On]/l3727 and [Om]AA4959, 5007 for run2 galaxies. At some 
distances from the center [Ora]/l4959 is not detected, in which 
case we assume [Om]i/l4959,5007=1.337[Oiii]^5007 (Oey & 
Kennicutt 1993). 

The abundance indicators used in this work are the follow- 
ing: 

1. R23: Based on this parameter {R23 = 
([<9//]i3727 + [OIII]AA4959, 5007])/ H/3), we use two 
calibrations: 

- M91: The theoretical calibration of McGaugh (1991), 
which takes into account the influence of the ionization 
parameter to determine the chemical abundance. 

- Z94: It is an analytical calibration obtained by Zaristky 
et al. (1994) as the average of three previous R23 cal- 
ibrations (Dopita & Evans 1986; McCall et al. 1985; 
Edmunds & Pagel 1984). It is only valid for high metal- 
licity objects with log(0/H)+12 > 8.35. 

2. N2: It uses ( m ^ m ) (Pettini & Pagel 2004) with the advan- 
tages of not being double valued with oxygen abundance, 
and that lines are very near each other, minimizing the de- 
pendence on the reddening correction. 

3. 03N2: It depends on ( ^"^m/Ha ) ( Pettini & p agel 2004), 
with the advantage of being a monotonic function of the 
oxygen abundance. However, in this case, the pairs of lines 
([O m]^5007, H/3) and ([N n]^6583, Ha) are detected in sep- 
arated spectra, which can introduce errors caused by a miss- 
match of the continuum fluxes in the red and blue (runl) 
ranges. It will be discussed at the end of this Section. 

4. 02Ne3: Perez-Montero et al. (2007) discussed the relation 
between [Ne iii]/13869 and [On], and define the parameter 
02Ne3 = mn^±l^Mimm _ Taking into account 
this parameter, the abundance is calculated with the expres- 
sion of M91. This line is more intense in AGNs due that, in 
these objects, Ne +2 is less suppressed than in star formation 
regions, as explained in Rola et al. (1997). 

5. SB1: It is calibrated for active galaxies (Storchi-Bergmann 
et al. 1998) and has two different calibrations. It is valid for 
Seyfert but not so accurate for LINERs. We use the calibra- 
tion involving the [Om], [Nn], Ha and H/3 lines because of 
our observed spectral range. It depends on the gas density, 
being valid only for 100<N< 10000 (with N the gas density 
in cm" 3 ) and 8.4 < 12 + log(0/H) < 9.4. 

6. KD02: Method proposed by Kewley & Dopita (2002) which 
consists on an iterative scheme involving the ionization pa- 
rameter and ^23 • 

7. P05: Pilyugin & Thuan (2005) improved the calibration of 
Pilyugin (2000, 2001) with a more extended (in metallicity) 
set of H 11 regions. The lines needed to measure this indicator 
are [On]AA 3727,3729, [Oni]^ 4959,5007 and H/3. 



R23, 02Ne3, KD03, P05 methods are applied to galaxies ob- 
served in run2 and N2, 03N2, SB1 to those observed in runl. 

Several of the indices described above show two branches 
in their calibration, i.e., there are two possible metallicity val- 
ues for a given index. This makes necessary to have further cri- 
teria to distinguish between the upper and the lower branches. 
We adopt the high-metallicity values on the basis of (a) the 
empirical diagrams of Nagao et al. (2006) (though not for ac- 
tive galaxies): With a sample of 50000 spectra they have pro- 
posed a method, based in the monotonic metallicity dependence 
of [Om]A 5007/[On]^ 3727, to break the R 23 degeneracy. We 
have made use of their diagrams to determine the branch, (b) 
the Kobulnicky et al. (1999) conclusion about galaxies in the lo- 
cal Universe (i.e. objects more luminous than Mb ~ -18 have 
metallicities higher than 8.3, and the whole sample lies within 
this range), (c) morphologically, all the galaxies are of early type 
and Oey & Kennicutt (1993) found that the average metallicity 
value for H11 regions in Sa and Sb galaxies is 8.97, clearly in 
the upper region, and (d) in NGC 1530 we obtain in this upper 
branch values similar to those found by Marquez et al. (2002) 
using the [Nul/Ha ratio. An special case is that of NGC 1358: 
the 7?23 values calculated for this galaxy are in the limit of appli- 
cability of the methods based on this parameter. But with data of 
[Nn] and Ha we have calculated abundance values supporting 
those obtained with R23 . 

The errors in the oxygen abundances have been estimated 
from error propagation on the uncertainties in the relevant line 
fluxes and range between ~0.01 and 0.10 dex (or 10-20% in 
[O/H]). This is lower than the intrinsic scatter associated to em- 
pirical calibrations of strong line methods (e.g. 0.2 dex for R23, 
Athey & Bregman 2009). We have carried out several tests to 
determine the uncertainties introduced in the different parts of 
the analysis. For that, abundances are estimated in two possible 
ways for each test of this analysis. In what follows, we describe 
the maximum difference found and the mean value of these dif- 
ferences in a galaxy, chosen as a characteristic example: 

(a) As discussed in Sect. 3, in most of the spectra we have used 
stellar templates to substract the contribution from the underly- 
ing stellar population, being a second possibility to correct the 
line intensities with an iterative method. The differences found 
in gaseous abundance between both methods are lower than 0.1- 
0.2 dex. As an example, for NGC 1433 the mean values of these 
differences are 0.04 + 0.06 dex for N2, 0.05 + 0.03 dex for 03N2 
and 0.04 ± 0.02 dex for SB1 

(b) Considering either an equivalent width due to Balmer absorp- 
tion of 2A to correct the stellar absorption, a typical value used 
in the literature (e.g. Zaritsky et al. 1994), or using the iterative 
method to estimate it, introduces differences lower than 0.1-0.2 
dex. For example, in NGC 1530 the mean values of these differ- 
ences are 0.04 + 0.04 dex (R 2i ), 0.05 + 0.04 dex (Z94), 0.03 + 
0.02 dex (KD03) and 0.13 ± 0.09 dex (P05) 

(c) The error introduced if we consider an l(Ha)/l(H/3) value of 
either 2.85 or 3.1 is as low as 0.05 dex. Annibaldi et al. 2010 
also showed that, in any case, considering HafH/3 — 2.85 does 
not significantly change the results. For NGC 3081, mean differ- 
ences are 0.0004 + 0.0008 dex (N2), 0.003 + 0.004 dex (03N2) 
and 0.08 + 0.09 dex (SB1) 

We next show that differences among the different metallicity 
calibrations are larger than these errors. 

We compare the results from the different methods in Fig. [8] 
taking N2 as reference for runl galaxies, and 7?23(M91) formn2 
galaxies. 02Ne3 can only be calculated for three galaxies; its 
behavior is the same as the /?23 (M91) but with a larger disper- 
sion. R 2i (Z94) is valid for 12 + log (O/H) > 8.35 (Zaritsky et 
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al. 1994). For 12 + log (O/H) values between 8.35 and 8.7, R 2i , 
with the calibrations of Z94 and M91, provide similar metallic- 
ity values, but from 8.7 upwards, the Z94 value becomes higher 
than the M91 one. In KD03 the behavior is the same as in R23 
(M91), but KD03 gives values slightly above those of R23 (M91) 
from 8.2 upwards. P05 values are systematically 0.3-0.5 dex 
lower than those of R 2 i (M91) with the exception of NGC 1358. 
Bresolin et al. (2009) showed that, for the external parts of the 
disk of M83, the metallicity with [Nii]/[Oii] is closer to the ac- 
tual values, being 0.4 dex lower than those obtained with R^i- 
Therefore, it might be expected that the metallicities using the 
P05 method are more realistic or, at least, closer to the direct 
metallicity determinations. Kennicutt et al. (2003) found system- 
atic differences up to 0.5 dex between T e and /?23 abundances. 
Pilyugin & Thuan (2005) noted that the determination of oxy- 
gen abundances using P05 agrees with T e abundance to within 
0. 1 dex, also for high-metallicity ones, which is our case. 

Forrunl galaxies, the largest differences between 03N2 and 
the N2 index are 0.3-0.4 dex, except for NGC 3081, in which 
case the difference is around 0.5 dex. Denicolo et al. (2002) 
emphasized the importance of the N2 method, showing that the 
03N2 calibration did not improve the results. SB1 is a method 
calibrated for AGN; though eleven of these galaxies have been 
classified as active, this method has only been used in the cases 
that meet the above conditions of density and metallicity, and 
only to runl galaxies because of the lines involved. 

It is interesting to compare the radial dependence of the val- 
ues obtained with the different methods, to study the effect of 
nuclear activity or the differences due to physical differencies 
in the nuclear and disk Hn regions properties. Hn characteris- 
tics can influence the abundance trends obtained using different 
methods in different ways (Kennicutt et al. 1989). The system- 
atic differences will be considered in the comparison with stellar 
metallicities. Figure [9] shows the residual values of the different 
abundance calculation methods for all the galaxies with respect 
to /?23 (M91) forrun2 galaxies and N2 forrunl galaxies. 

As already noted, P05 abundances are always below R23 
(M91), the difference being around 0.5 dex, without changes at 
different radii (with the expection of the special case NGC 1358 
as we have mentioned). For KD03, the largest differences are 
found in the central part of the Seyfert galaxy NGC 1358, being 
much smaller beyond 0.08 R/R25, the activity zone in Fig. [7] But 
this is not the general behavior for Seyferts (see NGC 2273), 
for which we find similar values to those calculated with R23 
calibrations by M91 and Z94. We find a larger dispersion for 
02Ne3 than for the other methods, without a systematic trend 
with radius, being the largest deviation for galaxies with less 
bins (e.g. NGC 4314). Between 03N2 and N2 the differences are 
small in most galaxies, less than 0.2-0.3 dex. But in two of them, 
NGC 3081 and NGC 2217, these differences are higher (0.4-0.5 
dex), and with the peculiarity that become smaller (again 0. 1-0.2 
dex) for larger galactocentric distances. 

As above mentioned, 11 galaxies have been classified as 
Seyferts or LINERs. All methods used here to derive metallicity, 
with the exception of SB1, have been derived for star-forming 
galaxies. We can only use the specific calibration method SB1 
in three out of our twenty galaxies: NGC 1433, NGC 2217 and 
NGC 3081, even only one of these, i.e. NGC 3081, is clearly a 
Sy2. Annibali et al. (2010) compared oxygen abundances esti- 
mated from a photoionization model based on R21 (Kobulnicky 
et al. 1999) with SB1. The values obtained from the SB1 method 
are on average ~0.04 dex larger than those using the R23 cal- 
ibration, but the differences between both metallicity values in 
individual regions can be as large as 0.3 dex. Keeping in mind 



this result we have opted to show the values obtained with all 
the methods for comparison. This can be specially useful in the 
case of LINERs, as for them the mechanism of ionization is not 
known. Therefore we have calculated the abundance values with 
SB1 when possible. For the active galaxies for which this method 
cannot be used we only give metallicities obtained from meth- 
ods calibrated for star forming regions, and we warn again the 
reader about the validity of these values in the AGN dominated 
area. 

4.2.2. Gas oxygen abundance in barred galaxies 

Figure [10] shows the oxygen abundances obtained with ^23 
(M91) or 03N2 versus radii for all galaxies. Linear fits, per- 
formed in the bar region and in the bulge region taking into ac- 
count the errors derived from the line fluxes, are also plotted for 
the different methods. 

Tables [3] and [5] show, respectively for run2 and runl, the 
value of oxygen abundances considering an aperture of 1.2 arc- 
sec (averaging abundances inside 1.2 arcsec), its central value 
from the extrapolation down to R=0 of the linear fit (only for 
7?23(M91) and 03N2), and the central value from the extrapola- 
tion down to R=0 excluding the influence region of the AGN. 
Tables [4] and [5] show gradients for all used methods (in units 
of dex/arcsec) taking into account all the nebular emission in 
the bar region and their errors. Galaxies classified as Seyferts or 
LINERs are marked with a dagger. For the galaxies for which 
we have enough points, we have two rows in these tables: The 
first one shows values obtained in the bar region, and the sec- 
ond one those obtained excluding the AGN influence region de- 
termined in the previous section. These influence regions have 
been determined from Fig. Qexcept for NGC 2681, NGC 2859 
and NGC 4394 for which we have considered the sizes reported 
by Veron-Cetty & Veron (2006) and references therein. The nu- 
merical values are shown in Tables|3]and|5] 

The seeing size could introduce errors in estimating the gra- 
dients. To study this effect we have produced artificial spec- 
tra at each pixel position. We have then convolved them with 
Gaussians of different width, covering the seeing values of these 
observations. This was made using the IRAF GAUSS task. We 
then generated three spectra. The largest differences in abun- 
dances calculated for these three spectra are lower than 0.2, im- 
plying differences in the gradient of ~ 5 %. 

We analyze the results using the indices 7?23(M91) and 
03N2. For run2 galaxies differences are nearly constant at all 
radii, in such a way that gradients maintain their sign in nearly 
all methods. In runl galaxies, there are changes with galactocen- 
tric distance as above mentioned. In particular, the three galaxies 
for which we have calculated SB1 coincide with those with the 
most different values of abundances in the central region, being 
smaller for 03N2 than for N2. This difference does not keep at 
larger radii, so the gradient sign is opposite for both methods. 
The oxygen abundances in a 1.2 arcsec aperture are in the range 
8.4-9.1 dex. 

If we consider all the detected emission in the bar region at 
a 3-cr significance, for the star formation galaxies, we detect a 
positive gradient for three galaxies (NGC 1 169, NGC 2523 and 
NGC 2935) and a negative gradient is observed for two galaxies 
(NGC 1832 and NGC 2962). For LINERs, we find a positive gra- 
dient for three galaxies (NGC 1433, NGC 4314 and NGC 4394) 
and none of them has a negative gradient. Considering the re- 
gion without influence of the possible AGN in these galaxies, 
the sign of the gradients does not vary. In the case of Seyferts 
we find negative gradients in four galaxies when we consider 
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Fig. 8. Relation between different calibrations for obtaining oxygen abundances. In (a), (b), (c) and (d) we compare the methods 
used for run2 galaxies with the /?23 calibration by McGaugh (M91). In (e) and (f) the comparison is for runl galaxies, between N2 
and 03N2 or SB1. The solid line in each plot shows the agreement between the different methods. 



the bar region, and excluding the AGN influence region, three of 
them maintain the gradient sign. Only one of the Seyfert galax- 
ies shows a positive gradient, which is not affected by the AGN 
region. 

In most of the disks of spiral galaxies, both barred or un- 
barred, negative metallicity gradients have been obtained (e.g. 
Bresolin et al. 2009, Garnett & Shields 1987, etc.). This gradient 



seems to be lower for barred galaxies than for unbarred ones, and 
higher for late type galaxies than for early type (Hidalgo-Gamez 
et al. 201 1, Vila-Costas & Edmunds 1992, Zaritsky et al. 1994). 
All our galaxies are barred and early type; so, based on previ- 
ous studies, the gradient should be shallow. We find gradients 
in the range [-3.7,1.4] (normalized to R25) with the exception of 
NGC 1169, which has a very high dispersion in the abundance 
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Fig. 9. Differences between the oxygen abundances obtained with the different methods versus the galactocentric radius for each 
galaxy. R25 has been taken from HyperLeda. Residuals are with respect to R23 (M91) for run2 galaxies and to N2 for runl galaxies. 
Crosses are oxygen abundances obtained from the 02Ne3 calibration, triangles from R23 (Z94), circles from KD03, squares from 
P05, diamons from 03N2 and asteriks from SB1. 
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Fig. 10. 12+log(0/H) obtained with 03N2 for run 1, and with R23 (M91) forrun2, versus R/R25 in each galaxy. Dashed lines indicate 
the fit for all measured values, and solid lines show the best fits within the bulge region. For runl galaxies, black lines indicate the 
N2, cyan lines the 03N2 fit and blue line the SB1 fit. For run2 galaxies black lines show the R21 (M91) fit, blue lines the Z94 fit, 
green lines the KD03 fit and red lines the P05 fit. 



Table 3. Central abundances, for run2 galaxies. 



Galaxy 


12+log(0/H) (1) 


12+log(0/H) (2) 


12+log(0/H)' 3) 


R (4) 
K 25 


AGN (5) 




dex 


dex 


dex 


arc sec 


arc sec 


NGC 1169 


8.47 ±0.10 


8.47 ± 0.01 




99.34 




NGC 1358t 


8.78 ± 0.07 


8.85 ± 0.04 


8.71 ±0.12 


64.14 


3.2 


NGC 1530 


8.90 ± 0.02 


9.071 ±0.006 




54.59 




NGC 1832 


9.09 ± 0.01 


9.075 ± 0.003 




73.64 




NGC 2273f 


8.57 ± 0.05 


8.809 ± 0.008 


8.828 ± 0.009 


68.72 


0.4 


NGC 2523 


8.48 ± 0.05 


8.50 ±0.10 




84.55 




NGC 268 If 


8.92 ± 0.02 


9.007 ± 0.007 


9.12 ±0.02 


119.43 


3* 


NGC 2859f 


8.66 ± 0.03 


8.872 ± 0.007 


9.05 ± 0.01 


94.87 


3* 


NGC 2962 


8.58 ± 0.03 


8.72 ± 0.01 




68.72 




NGC4314f 


8.81 ± 0.02 


8.900 ± 0.004 


8.965 ± 0.008 


111.46 


3.3 


NGC 4394f 


8.69 ± 0.03 


8.708 ± 0.009 


8.67 ± 0.02 


104.02 


3* 



tActive galaxies. 

'''Central oxygen abundance in the gas, estimated with an aperture of 1.2 arcsec, and calculated with the R 2 3 (M91) method. 

(2) Central oxygen abundance in the gas, as a result of the fit, calculated with the R23 (M91) method. 

(3) Central oxygen abundance in the gas, as a result of the fit outside the region of AGN influence, calculated with the ^23 (M91) method. 

(4) R25 is from the Hyperleda galaxy catalogue. 

(5) AGN influence region, in arcsec, obtained from Fig.|7]except those marked with *, which values are obtained from Veron-Cetty & Veron (2006) 
and references therein. 



values. In Sect. 5 we compare these values with those found in 
the literature. 

With respect to the bulge region, in Tables |6]and |7]are shown 
the obtained gradients for all used method, and the central oxy- 
gen abundance (only with R23 (M91 ) or 03N2 method). In these 



tables we have marked the active galaxies but the gradients are 
fitted including the AGN regions. NGC 3081 is the only Seyfert 
galaxy for which we have been able to calculate abundances with 
the SB1 method. For two galaxies (NGC 1169 and NGC 5101) 
all the measured nebular emission in the bar position is in the 
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Fig. 10. Continued 



Table 4. Oxygen abundance gradients with all used methods and the central abundance, for run2 galaxies. 



Galaxy 


grad-fi 23 (M91) 


grad-02/Vei 


grad-^ 23 (Z94) 


gmd-KD03 


grad- P05 




dex/arcsec 


dex/arcsec 


dex/arcsec 


dex/arcsec 


dex/arcsec 



NGC 


1169 


0.080 ± 0.007 




0.108 ±0.015 


0.086 ± 0.012 


0.055 ±0.013 


NGC 


1358t 


-0.058 ± 0.017 


0.04 ± 0.03 


0.11 ±0.01 


-0.05 ± 0.01 


-0.12 ±0.05 






-0.03 ± 0.02 


0.06 ± 0.08 


0.05 ± 0.03 


-0.03 ± 0.02 


-0.08 ±0.10 


NGC 


1530 


-0.0006± 0.0002 




0.0014 ± 0.0002 


-0.0014 ± 0.0003 


0.0069 ± 0.0002 


NGC 


1832 


-0.0040 ± 0.0002 




-0.0036 ± 0.0003 


-0.0005 ± 0.0002 


-0.0102 ± 0.0002 


NGC 


2273f 


-0.0072 ± 0.0013 


-0.01 ± 0.02 


-0.001 ± 0.002 


-0.0055 ± 0.0015 


-0.012 ± 0.002 






-0.009 ± 0.001 


-0.07 ± 0.02 


-0.004 ± 0.002 


-0.006 ± 0.001 


-0.014 ± 0.002 


NGC 


2523 


0.0093 ± 0.0004 




0.0133 ± 0.0008 


0.0073 ± 0.0005 


0.0045 ± 0.0007 


NGC 


268 It 


-0.020 ± 0.002 




-0.014 ± 0.004 


-0.012 ± 0.001 


-0.029 ± 0.003 






-0.046 ± 0.005 




-0.051 ± 0.009 


0.010 ± 0.004 


-0.074 ± 0.008 


NGC 


2859t 


-0.0078 ± 0.0013 




-0.001 ± 0.003 


0.000 ± 0.002 


-0.012 ± 0.003 






-0.027 ± 0.002 




-0.028 ± 0.004 


-0.020 ± 0.003 


-0.038 ± 0.003 


NGC 


2962 


-0.017 ± 0.002 




-0.026 ± 0.004 


-0.016 ± 0.004 


-0.018 ±0.003 


NGC 


4314t 


0.0123 ± 0.0007 


0.111 +0.014 


0.0165 ± 0.0014 


0.0066 ± 0.0007 


0.0136 ±0.0012 






0.005 ± 0.001 




0.0016 ± 0.0019 


0.0015 ± 0.0012 


0.0046 ± 0.0016 


NGC 


4394f 


0.0073 ± 0.0003 




0.0097 ± 0.0005 


0.0051 ± 0.0003 


0.0068 ± 0.0004 






0.0080 ± 0.0008 




0.0104 ± 0.0014 


0.006 ± 0.001 


0.0084 ± 0.0012 



tActive galaxies (Seyferts and LINERs): The first row shows values calculated with all points, and the second one those calculated in the region 
outside the AGN influence. 



bulge zone where the scatter is larger than in the bar region and, 4.3. Comparison of the nebular gas abundances with the 
therefore, errors are also larger. We find the same sign for the distribution of stellar metallicities 

slopes in most of the star-forming galaxies. They are in general 
steeper than those considering all the bar region. 

While measurements of the ionised gas abundances give infor- 
mation about the current composition of the ISM, the stellar 
metallicities are directly linked to the star formation history and, 
thus, a comparison between the metal abundances of both com- 
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Table 5. Oxygen abundance gradients with all the used methods, forrunl galaxies. Units are dex/arcsec. 



Galaxy 


grad-yV2 
dex/arcsec 


gmd-03N2 
dex/arcsec 


grad-Sfii 
dex/arcsec 


124-log(0/H) (1) 
dex 


124-log(0/H) (2) 
dex 


arcsec 


AGN 
arcsec 


NGC 1433t 


-0.0022 ± 0.0002 


0.0012 4 


: 0.0002 


0.01 ± 


0.01 


8.68 ± 


0.03 


8.724 = 


t 0.008 


184.98 






-0.0014 ± 0.0002 


0.0009 i 


: 0.0002 










8.77 = 


tO.Ol 




3.7 


NGC 2217f 


-0.0054 ± 0.0009 


0.0021 i 


: 0.0008 


-0.023 d 


: 0.007 


8.674 ± 


0.03 


8.699 = 


t 0.009 


140.32 






0.0004 ± 0.0025 


-0.0002 : 


± 0.0009 










8.76 = 


t0.03 




16.8 


NGC 2665 


-0.0041 ± 0.0006 


-0.0010 : 


± 0.0006 






8.718 ± 


0.03 


8.739 = 


t 0.009 


59.86 




NGC 2935 


-0.0034 ± 0.0003 


0.0019 4 


: 0.0004 






8.649 ± 


0.02 


8.714 = 


t 0.007 


125.06 




NGC 308 It 


-0.0023 ± 0.0010 


0.0142 4 


: 0.0009 


0.024 ± 


0.012 


8.369 ± 


0.02 


8.332 = 


t 0.009 


75.36 






-0.0025 ± 0.0021 


0.0187 4 


: 0.0018 










8.22 = 


t0.04 




7.5 


NGC 4643 f 


0.004 ± 0.003 


-0.003 : 


± 0.002 






8.751 ± 


0.06 


8.756 = 


t 0.019 


67.16 


13.2 


NGC5101f 


-0.028 ± 0.008 


-0.004 : 


± 0.008 






8.784 ± 


0.03 


8.780 = 


t 0.019 


176.66 


3.5 



'''Central oxygen abundance in the gas, estimated with an aperture of 1.2 arcsec, and calculated with the 03N2 method. 

* 2 'Central oxygen abundance in the gas, as a result of the fit, calculated with the 03N2 method. 

* 3 'R25 is from the Hyperleda galaxy catalogue. 

< 5 'AGN influence region, in arcsec, obtained from Fig. [7] 

tSame as Tableg] 



Table 6. Metallicity gradients in the bulge regions with all the used methods forrun2 galaxies, including AGN regions. 



Galaxy 


grad-^3 (M91) 


12+log(0/H)"' 


grad-02AW 


grad-^sf Z94) 


giad-KD03 


grad-P05 


order' 2 ' 




dex/arcsec 


dex 


dex/arcsec 


dex/arcsec 


dex/arcsec 


dex/arcsec 




NGC 1169 


0.080 ± 0.007 


8.47 ± 0.01 




0.108 ±0.015 


0.086 ± 0.012 


0.055 ± 0.013 


1 


NGC 1358f 


-0.08 ± 0.02 


8.88 ± 0.05 


0.05 ± 0.04 


0.16 ±0.03 


-0.069 ± 0.019 


-0.15 ±0.10 


2 


NGC 1530 


0.008 ± 0.006 


9.03 ± 0.02 




0.016 ±0.013 


0.015 ± 0.009 


0.023 ± 0.006 


3 


NGC 1832 


-0.038± 0.003 


9.113 ±0.007 




-0.053± 0.005 


-0.020 ± 0.003 


-0.053 ± 0.004 


4 


NGC 2273f 


-0.095± 0.007 


8.892 ± 0.014 


-0.01 ± 0.02 


-0.117+ 0.012 


-0.078 ± 0.007 


-0.08 ± 0.01 


5 


NGC 2523 


0.054 ± 0.010 


8.58 ± 0.02 




0.082 ± 0.017 


0.058 ±0.011 


0.053 ± 0.014 


6 


NGC 268 If 


-0.01 1± 0.003 


8.987 ± 0.007 




-0.005± 0.005 


-0.011 ±0.001 


-0.018 ± 0.004 


7 


NGC 2859t 


0.037 ± 0.003 


8.762 ± 0.009 




0.046 ± 0.006 


0.024 ± 0.003 


0.036 ± 0.006 


8 


NGC 2962 


0.109 ±0.006 


8.48 ± 0.01 




0.111 ±0.012 


0.104 ±0.009 


0.178 ±0.010 


9 


NGC4314f 


0.0189 ±0.0012 


8.871 ± 0.005 


-0.11 ±0.16 


0.034 ± 0.002 


0.0192 ± 0.0015 


0.0106 ± 0.0019 


10 


NGC 4394f 


0.035 ± 0.010 


8.671 ±0.017 




0.044 ± 0.017 


0.042 ±0.011 


0.053 ±0.015 


11 



'"Central oxygen abundance in the gas, as a result of the fit in the bulge region, calculated with the R03 (M91) method. 
(2) SeeFigs.[ni& [13] 
tSame as TableU 



Table 7. Metallicity gradients in the bulge regions, with all the used methods, forrunl galaxies, including AGN regions. 



Galaxy 


grad-JV2 
dex/arcsec 


12+log(0/H)'" 
dex 


grad-OJA^ 
dex/arcsec 


gmd-SBl 
dex/arcsec 


order' 2 ' 


NGC 1433t 


-0.001 : 


± 0.006 


8.665 


± 


0.017 


0.024 ± 0.006 


0.01 4 


:0.01 


12 


NGC 2217f 


-0.045 : 


± 0.006 


8.698 


± 


0.018 


0.003 ± 0.005 


-0.023 : 


±0.012 


13 


NGC 2665 


0.008 i 


: 0.007 


8.709 


± 


0.018 


0.012 ± 0.006 






14 


NGC 2935 


0.005 4 


: 0.004 


8.653 


± 


0.016 


0.015 ± 0.004 






15 


NGC 308 It 


0.002 4 


: 0.005 


8.376 


± 


0.015 


-0.003± 0.005 


0.022 ± 0.015 


16 


NGC 4643 1 


0.009 =t 


: 0.004 


8.73 


± 


0.02 


0.004 ± 0.004 






17 


NGC5101f 


-0.028 : 


± 0.008 


8.780 


± 


0.019 


-0.004± 0.008 






18 



'"Central oxygen abundance in the gas, as a result of the fit in the bulge region, calculated with the 03N2 method. 

< 2) SeeFigs.[[2]& E2] 
tSame as Tableg] 



ponents can help us to determine the origin of the gas involved 
in the present star formation. 

In previous works (Perez et al. 2007, 2009, 2011) we de- 
rived the radial distribution of stellar ages and metallicities along 
the bar and the bulges of the galaxies presented in this work. 
Line strength indices were measured and used to derive age and 
metallicity gradients in the bulge and bar region by comparing 
with stellar population models. 



Figure QT| shows the comparison of the stellar and ionized 
gas abundances along the radius for all the galaxies. For the sake 
of clarity the radius is on a logarithmic scale and R=0 is shifted 
to log (R)=-0.5. Forrun2 galaxies the Rri(M9\) calibration has 
been chosen, covering the error bars the values obtained with all 
methods with the exception of 02Ne3. The values obtained with 
all methods have been represented in the case of runl galax- 
ies. For the comparison, we take as reference the solar value 
12+log(0/H) =8.69 (Allende-Prieto et al. 2001) and [Fe/H] for 
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both components. It has been also shown the AGN character and 
we have marked with a vertical line the influence region of these 
active nuclei. In this way, we warn again (Sect. 4.2.1) the reader 
about comparing the gas abundances with the stellar metallici- 
ties in the regions dominated by an AGN. 



For the following mn2 galaxies the nebular metallicities 
lie clearly below the stellar metallicity values: NGC 1169, 
NGC 2962 (in agreement with Marino et al. 2011) and 
NGC 4394. Only NGC 4394 is classified as LINER. This lower 
metallicity gas (w.r.t. the stellar metallicity) is located in the 
bulge region, similar values for the metallicity of both compo- 
nents are found elsewhere. For runl galaxies nebular metallici- 
ties obtained with 03N2 are smaller than stellar metallicities in 
NGC 2217, NGC 3081 and NGC 4643. But with the values ob- 
tained with N2 they are very similar (NGC 308 1 and NGC 4643) 
or even higher (NGC 2217). These three galaxies are classified 
as AGNs. If we consider SB1, nebular metallicity is of the same 
order than stellar metallicity for NGC 3081 (Sy) and NGC 2217 
(LINER). 

For some of the other galaxies the distribution of the gas and 
the stellar abundances closely follow each other and for 3 galax- 
ies the gas abundances lie clearly (>0.2 dex) above the stellar 
metallicities. 

Not many works have addressed a direct comparison be- 
tween the metallicities obtained from the stars and the ionized 
gas. Storchi-Bergmann et al. (1994) concluded that both abun- 
dance measurements are well correlated for star-forming galax- 
ies. However, recent work (Marino et al. 2011) suggests, as de- 
rived from HI and GALEX UV imaging, the presence of current 
star formation activity in early type galaxies with emission lines. 
This sort of galaxy rejuvenation could be due to external gas re- 
fueling. This accretion would suggest differences in the metal- 
licities of the old stellar component and the newly accreted gas, 
responsible for the present star formation. Furthermore, Annibali 
et al. (2010), found in a comparison between nebular and stellar 
metallicities of E-S0 galaxies with nebular emission that the gas 
metallicity tends to be lower that the stellar metallicity, being the 
effect most severe for the largest stellar metallicities. 

In Figs. [T2] & [13] we compare the gradients obtained 
for the oxygen abundances along bars and the bulge regions 
with those estimated in Paper II. They are normalized to R25 
((12+log(0/H))/(R/R 25 )). In Paper I we found that the stellar 
metallicity gradients along the bar show a large variety. The dis- 
tribution of the gas nebular abundances in the bar region shows 
a mild gradient (flatter than for the stellar component for star- 
formation galaxies except for NGC 1169), and we do not see 
any clear relation between the stellar and nebular gas metallic- 
ity gradients. With respect to the bulge region, gradients for the 
stellar component in star-formation galaxies are flatter than the 
gaseous gradients for NGC 1169 and NGC 2523. In the mea- 
surable cases, the oxygen abundance gradients are, in general, 
steeper for bulges than for bars, and they have the same sign. 
However, this is not the case for the stellar abundance gradient, 
where the sign between the bulge and the bar gradient changes 
in 10 out of 16 galaxies. 

5. Discussion 

Most of the emission is centrally concentrated in the bulge re- 
gion. This has been noted before and seems to be independent of 
the presence of a bar (Davis et al. 201 1). In 5 galaxies, emission 
lines are present along the whole extension of the bar, and 6 out 
of the 18 galaxies present some star formation at the bar ends as 
well as in the bulge region. Previous works have claimed that the 
distribution of star formation in a bar is correlated with the bar 
age; younger bars (< lGyr) show star formation along the entire 
bar while older bar stages show star formation at the bar-ends 
and nuclear region (Phillips 1993; Martin & Friedli 1997; Martin 
& Roy 1995; Friedli & Benz 1995; Verley et al. 2007). However, 
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Fig. 12. Comparation between metallicity gradients obtained for 
gaseous and stellar component in the bar region of our sample 
galaxies. They are normalized to R25 (i.e. the gradients have been 
obtained from linear fits to the 12+log(0/H) versus R/R25 repre- 
sentation for each galaxy). Circles are those estimated with the 
7?23(M91) or N2 methods, squares are the mean value of gradi- 
ents calculated with all methods (the error bar covers the range 
of these gradients), and stars are the gradient [Z/H] for the stellar 
component (Paper II). Order in abcisas is that of Tables [6]and [7] 
Filled symbols are for star-formation galaxies. 



some of the galaxies for which nebular gas is found along the bar 
seem to host bars that are old (see Paper I). Therefore, the star 
formation distribution in these bars is possibly more related to 
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Fig. 13. Comparation between metallicity gradients obtained for 
gaseous and stellar component in the bulge region of our sam- 
ple galaxies. They are normalized to R25 (i.e. the gradients have 
been obtained from linear fits to the 12+log(0/H) versus R/R25 
representation for each galaxy). Circles are those estimated with 
the 7?23(M91) or N2 methods, squares are the mean value of gra- 
dients calculated with all methods (the error bar covers the range 
of these gradients), and stars are the gradient [Z/H] for the stellar 
component (Paper II). Order in abcisas is that of Tables [6]and [7] 
Filled symbols are for star-formation galaxies. 



the reservoirs of gas that could replenish the bar than to the ac- 
tual time of bar formation. 

We do not find neither a relation between the difference in 
abundance values and the type of the galaxy nor a relation with 
the galactocentric distance in which the lines are measured (see 
Sect. 4.2.1 and 4.2.2). 

We have considered the possible systematic effect of the im- 
pact of an AGN nuclei in the determination of abundances. Its 
effect is difficult to quantify, but very important when compar- 
ing with the stellar metallicities in the inner regions of these 
galaxies and deriving the abundance gradients. The central re- 
gions will be a mixture of different components with different 
physical properties (Lindblad & Fathi 2011) and will introduce 
difficulties in determining the main ionising mechanisms. Most 
of the central abundances in the works above mentioned, have 
not considered the influence of an active nuclei. Their abun- 
dance estimations have been carried out in H 11 regions or/and in 
star formation galaxies, so the models for calibration are based 
on stellar photoionization alone. Storchi-Bergmann et al. (1998) 
proposed two calibrations for active galaxies: the first one de- 
pending on [NrfJ/Ho' and [Om]/HyS (SB1), and the second, a lin- 
ear combination of log [Nn]/Ha' and [On]/[Om] (SB2). They 
concluded that, when possible, the best is to use the mean value 
of both. Based on the extrapolation of metallicities towards the 
nuclear regions, they obtained that this calibration is good for 
Seyferts, but there is no clear conclusion for LINERs. We have 
used this method, SB1, only for three galaxies of runl. We can- 
not estimate SB2 for them because we cannot measure [0 11] . For 
these galaxies, 2 LINERs and a Sy, we compare abundances cal- 
culated with SB1 and N2 and 03N2, that would be the extreme 
cases: i.e. to consider AGN or photoionization produced by hot 
stars. For LINERs, the results of SB1 matches or is closer to that 
of 03N2 in the central regions, while for the Seyfert galaxy they 
are similar to those estimated with N2. In any case, the difference 
between oxygen abundance using N2 or SB1 is less than 0.2 dex, 
a value considered as the intrinsic uncertainty of the models used 
to calibrate R^i- Although the presence of an AGN does not seem 
to alter the results, we would like to warn the reader about the 
validity of the results in the regions influenced by active nuclei. 

We compare oxygen abundances with the stellar metallic- 
ities. We find central metallicities, 12+log(0/H), with a mean 
value within 1.2 arcsec, in the range of 8.4 to 9.1. These val- 
ues agree with those found in the inner parts of early-type spi- 
rals. Oey & Kennicutt (1993) studied Hn regions in 15 Sa to Sb 
galaxies, obtaining abundance values, with R^i, from 8.7 to 9.2. 
The larger values correspond to those regions located at lower 
galactocentric radii. They find that gas metallicities in these 
galaxies are systematically higher than those in Sc and later type 
galaxies. Zaritsky et al. (1994) discuss that this relation might 
reflect a dependence with galaxy mass and not with galaxy type 
after analysing H 11 regions in a sample of 39 galaxies covering 
morphological types from Sab to Sd. They found metallicities 
in the range 8.2-8.6, similar to the abundances of our work, al- 
though we cover morphological types from SBO to SBbc. Vila- 
Costas & Edmunds (1992) also find a correlation between cen- 
tral gas abundances and galaxy mass, analysing the gas metal- 
licity distribution of 32 discs of different morphological types. 
For the inner radii of Sab-Sb they find oxygen abundances in 
the range 8.7-9.5. A large statistical study using SDSS data by 
Tremonti et al. (2004) analyzes 53.000 star forming galaxies at 
z ~ 0. 1 . In this study, the authors address the fundamental role 
that the galaxy mass plays in its chemical evolution. They find 
a strong correlation between mass and nebular abundance for 
galaxies with masses between 10 8 5 and 10 M , with a flatten- 
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Fig. 14. Rotational velocity corrected for inclination versus cen- 
tral oxygen abundances (1.2 arcsec). Filled circles represent 
those calculated using the R23 method for run2 and 03N2 for 
runl. Squares are values obtained with N2 for run 1. 
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Fig. 15. Stellar velocity dispersion versus central oxygen abun- 
dances (1.2 arcsec). Filled circles represent those calculated us- 
ing the R23 method for mn2 and 03N2 for runl. Squares are 
values obtained with N2 for runl . 



ing of the relation for larger galaxy masses. Some recent studies 
investigating how mergers and bars can modify this correlation 
(Ellison et al. 2011) find a larger nebular abundance in barred 
galaxies than in unbarred ones of the same mass. Our galaxies 
lie on the flat part of the correlation, and therefore we do not ex- 
pect any strong correlation between the rotational velocity and 
the nebular abundance. 

However, we see in Fig. [14] that the galaxies with the low- 
est abundances are those with the largest rotational veloci- 
ties, although there is a large dispersion in the values, also 
finding galaxies with high rotational velocity and large ion- 
ized gas metallicities. However, there are no galaxies with low 
(12+log(0/H) < 8.68) nebular abundances and low rotational ve- 
locities. The only exception is NGC 3081 (Seyfert), but with the 
abundance value obtained using 03N2. This result might be re- 
flecting that it cannot be rejected some rejuvenation mechanism 
for some of the most massive galaxies in our sample. 

Following this result on the tendency of nebular abundances 
with disk rotational velocities, we search for nebular abundance 
trends with stellar velocity dispersion; however, we find no such 
trend. As for the stellar abundances, we find in paper II, as ex- 
pected from other studies, a correlation between stellar central 
abundances and central velocity dispersion (see Fig. Q3). 

Regarding the comparison of these central nebular abun- 
dances with central stellar metallicities, we find that for three 
of the galaxies the nebular abundances lie well bellow the stel- 
lar metallicities. The galaxies which show nebular abundances 
lower than the stellar metallicities are systematically among 
those galaxies with the largest rotational velocities of our sam- 
ple. 

The bar is thought to be an efficient mechanism to quickly 
redistribute material within the disk. External accretion via the 
bar potential is a plausible mechanism to bring new gas to the 
center and dilute the central gas metallicities. The strength of 
the bar determines the gas inflow rate (e.g. Athanassoula 1992, 
2000); however, no clear trend between the bar strength and neb- 
ular composition is found for our galaxies. The efficiency of this 
process might be limited to the availability of gas, and although 
no trend with HI mass is found we do not posses a complete in- 
formation about the total gas content of our galaxies. Other pro- 



cesses capable of explaining the existence of low metallicity gas 
in the center of these massive galaxies include cold flows to the 
center (Keres et al. 2005) and particular star formation histories 
(Yates et al. 2011). Yates et al. (2011) conclude, by examining 
semi-analytical models, that massive galaxies with low nebular 
abundances have undergone a gas-rich merger with a later shut- 
down of star formation. This low density metal-poor gas accre- 
tion is not sufficient to form stars efficiently and subsequently di- 
lutes the gas-phase metallicities. Further information of the total 
gas content, its distribution and kinematics would help to eluci- 
date the different external accretion mechanisms that might be at 
play in these galaxies. 

As for the nebular abundance gradients found in this work, 
numerical models of barred galaxies (Friedli et al. 1994) show 
that the initial slope of the abundance is only slightly modified 
in the bar region, while outside the bar corotation radius, the 
original slope gets flatten. We have estimated the gas oxygen 
abundances along the bars to study the efficiency of mixing by 
the bar inside corotation. Considering only star-forming galax- 
ies, we find positive gradients, with a 3-<x significance, for the 
oxygen abundance in the bar region in three galaxies and neg- 
ative gradients in two of the galaxies. The estimated gradient 
values range from -3.7 to 1.4, in units of dex/(R/R25) with the 
mentioned exception of NGC 1 169. In 8 of the observed galaxies 
the gradient is lower than 0.5 dex/(R/R25). In the bulge region, 
where most of the emission is found, we find steeper gradients. 

It is, however, not a simple task to interpret the origin of 
these gradients because not only stellar evolution, star formation 
and bar dynamics, already complicated processes, affect the gra- 
dients but also refueling of newly accreted gas could change the 
gradient. Naively, central accretion of gas would dilute the inner 
regions generating a positive gradient, while the effect of radially 
accreted gas could be not so evident without proper modelling. 
However, it would be possibly to reduce a pre-existing gradient 
because of the rapid mixing and reactivation of star formation in 
the bar region (Friedli et al. 1994, 1995). 

Most of the gas metallicity gradients that can be found in 
the literature include the whole disk; van Zee et al. (1998) found 
gradients in the range [-0.30, -1.52] dex/R25 for spiral galax- 
ies, and Rupke et al. (2010) a mean value of -0.57 + 0.05 con- 
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sidering 11 isolated spirals. Oey & Kennicutt (1993) concluded 
that barred galaxies have a shallower gradient than unbarred 
ones. But a large range of abundance gradients have been found, 
some of them being much more steeper than for unbarred spirals 
(Edmunds & Roy 1993, Considere et al. 2000). Vila-Costas & 
Edmunds, found gradients raging from -3.33 to -0.899 for non- 
barred galaxies, and from -0.156 to -0.398 for barred galaxies 
(all in dex/R25). Zaritsky et al. (1994) obtained a mean value 
of -0.23 dex/R25 considering five barred galaxies. In particular, 
for the Milky Way, Balser et al. (201 1) find a radial gradient of 
-0.0446 + 0.0049 dex/kpc. Advardsson (2002) found a value of 
-0.07 + 0.01 dex/kpc for galactocentric distances between 6 and 
18kpc. 

The discrepancy between the values found in the literature 
for the disk and our values, concerning only the bar region, could 
be due to the predicted change in the slope around the bar radius 
(Friedli et al. 1994). Furthermore, we find a slope change be- 
tween the bulge and the bar region, also predicted by Friedli et al. 
(1994). These changes in the abundance slope with radius, bee- 
ing steeper in the inner parts of the galaxies, have been analised 
by other authors. Zaritsky et al. (1994) found different gradient 
in the bar of NGC 3319 than outside of it. This question has 
been studied in later works, trying to detect a radial variation 
of the gradient. For instance, for the barred galaxy NGC 3359 
Martin & Roy (1995) found a radial variation of the gradient 
being, in turn, studied by Zahid & Bresolin (2011), who esti- 
mated a break at a characteristic radius. Balser et al. (201 1 ) con- 
cluded that there was no radial discontinuities for the Milky Way, 
but that the gradient varies in the azimuthal coordinate between 
-0.03 and -0.07 dex/kpc. 

Previous works have assessed the effect of the AGN pres- 
ence in the nebular abundance gradients. Annibali et al. (2010) 
showed that the nebular abundances, derived with the R23 
method, increase with radius, while they decrease when us- 
ing methods that take into account harder ionising sources (i.e 
AGN), such as the SB method, a combination of SB1 and SB2. 
We have checked whether our gradient calculation could be bi- 
ased in the same way. We can estimate SB 1 for three runl galax- 
ies, NGC 1433 and NGC 2217 and NGC 3081. For them, we 
have not [On], so we cannot calculate Rzs, but N2 and 03N2. 
We find that SB1 gradients are nearly flat for NGC 1433 and 
NGC 3081, and have the same sign than that from N2 for 
NGC 2217. 

There have been studies correlating the observed nebular 
metallicity gradients with physical properties of the galaxies. 
Zaritsky et al. (1994) suggested a correlation between the gradi- 
ent and the bar-type as previously noticed by Pagel & Edmunds 
(1981), Martin (1992) and Edmunds & Roy (1993). We have 
not found any correlation between the estimated gradients and 
the physical properties of galaxies: e.g.; bar-type, HI mass or 
central stellar velocity dispersion. However, we cover a short 
range of masses and morphological types. Further studies cover- 
ing a larger range of galaxy parameters and spatial range will be 
necessary to draw strong conclusions about the correlation be- 
tween nebular abundance gradients and the physical properties 
of galaxies. 

6. Summary and conclusions 

We have carried out a detailed analysis of the nebular abun- 
dances along the bar and in the bulge of a sample of 20 early- 
type galaxies to compare them with the stellar abundance distri- 
butions obtained in previous works using line-strength indices. 
We focus on several relations, in particular, the distribution of 



abundances and star formation and the properties of the galax- 
ies; the interpretation of the ionizing mechanisms and the rela- 
tive influence of the AGN; the relation between abundances and 
rotational velocities and velocity dispersion in the bulge. We 
have found that most of the emission in our sample of galax- 
ies are concentrated in the bulge region. We have used sev- 
eral methods to estimate nebular abundances. Considering dif- 
ferent methods for star-forming galaxies, does not change our 
results about gaseous metallicity gradients nor change the com- 
parison with the stellar metalicities. In three galaxies the esti- 
mated nebular abundances lie clearly below the stellar metal- 
licities. Furthermore, we see that the galaxies with the lowest 
abundances are those with the largest rotational velocities, and 
although there is a large dispersion in the values, this effect 
might be the result of a central rejuvenation mechanism in the 
most massive and late-type galaxies. The comparison between 
gaseous and stellar component shows that in most galaxies the 
slope for the nebular abundance distribution in the bulge region 
is shallower than that of the stellar metallicity. We can con- 
clude that the combination of observations of both gas and stellar 
metallicities is crucial to determine the star formation history in 
galaxies. The study of both phases is very useful to learn about 
the origin of the gas involved in the star formation; for instance, 
we clearly see galaxies for which external accreted gas is the fuel 
source for the present star formation. Due to the small parame- 
ter range covered by our sample we can not conclude anything 
about the role of bars in the mixing, but we observe a variety of 
gas and stellar metallicities distributions that indicate very differ- 
ent mixing and inflow histories from galaxy to galaxy. Therefore, 
it would be very interesting to complete this study with a wider 
sample, including later-type and unbarred galaxies. 
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Fig. 11. Oxygen gas abundance and stellar metallicities relative to solar as a function of logarithm of galactocentric radius in arcsec. 
R= is shifted to log/? = -0.5. Asterisks indicate metallicities relative to solar obtained from the stellar population (see text for 
details). Open squares indicate oxygen gas abundances using R23 for run2 galaxies, and error bars cover values obtained with all 
methods but 02Ne3. Fon runl galaxies, open squares indicate those obtained with 03N2, diamonds with N2 and triangles with SB1. 
AGN galaxies are marked, with the vertical lines indicating the influence region of the active nuclei (either as in Fig. [7] or 3", as 
indicated in references in Veron-Cetty & Veron (2006)). 
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Fig. 11. Continued. 
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Fig. 11. Continued. 
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